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We study the mobility of Dirac fermions in monolayer graphene on a GaAs substrate, restricted 
by the combined action of the extrinsic potential of piezoelectric surface acoustical phonons of GaAs 
(PA) and of the intrinsic deformation potential of acoustical eigen-phonons in graphene (DA). In the 
high temperature (T) regime the momentum relaxation rate exhibits the same linear dependence on 
T but different dependences on the carrier density n, corresponding to the mobility /i oc l/y/ri and 
1/n, respectively for the PA and DA scattering mechanisms. In the low T Bloch-Griineisen regime, 
the mobility shows the same square-root density dependence, [i oc •Jn, but different temperature 
dependences, fj, oc T -3 and T -4 , respectively for PA and DA phonon scattering. 



Graphene pQ due to its unique linear chiral electronic 
dispersion [2] exhibits novel transport properties [31 |3] 
and has great potential as a desirable material for future 
electronic and optical technologies [H [5] . Momentum re- 
laxation is a key phenomenon that governs transport of 
Dirac fermions in graphene [6] . It is of practical interest 
for developing high-speed electronics and in recent years 
has been extensively studied both theoretically [7HH] and 
experimentally [TUHE] ■ The scattering by defects [7J IT51 - 
[T5], impurities E2 QTHH] , and phonons HQ E3 [H- 
[26] have been investigated to determine and control the 
dominant mechanism that limits the carrier mobility in 
graphene. 

In device structures used so far, graphene is often de- 
posited on an oxidized silicon wafer (SiO^/Si), which due 
to various scattering mechanisms imposes constraints on 
its excellent transport properties observed in suspended 
graphene devices fTUl H7J ■ Recently, structures on other 
promising substrate materials such as h-BN [28l [29] and 
GaAs [3D1 [31] have been fabricated and studied with the 
intention for high-quality graphene electronics. Along 
with its superior surface quality and strong hydrophilicity 
preventing folding of large-scale graphene flakes, GaAs 
has a substantially larger dielectric constant in compar- 
ison with SiO*2 and h-BN and hence improved electrical 
screening. In such high purity GaAs structures, electron- 
phonon scattering can be a decisive factor in limiting the 
mobility of Dirac fermions and the piezoelectric GaAs 
substrate can serve as a powerful tool for studying elec- 
tronic properties of graphene by means of remote piezo- 
electric surface acoustical phonons. 

In the present work we study the temperature and 
density dependence of the carrier mobility in monolayer 
graphene at finite doping on a GaAs substrate. We cal- 
culate the mobility limited by scattering from the piezo- 
electric potential of remote surface acoustical phonons of 
the substrate (PA phonons) versus the deformation po- 
tential of acoustical eigen-phonons of the graphene lattice 
(DA phonons). In experiment the typical wavelength of 



phonons taking part in scattering events is much larger 
than the distance, d, of several angstroms between the 
graphene sheet and the GaAs substrate so that extrinsic 
interaction of remote PA phonons with Dirac fermions 
is quite strong and, as we shall see, can dominate the 
intrinsic interaction with the DA phonons of graphene. 

Crystal surfaces modify substantially bulk phonon 
modes and can change qualitatively the carrier relax- 
ation characteristics [32] ■ In crystals with lack of a cen- 
ter of symmetry such as in GaAs, the displacement field 
of Rayleigh waves (a combination of longitudinal and 
transversal oscillations) propagating on a crystal free sur- 
face [33) induces a piezoelectric polarization of the lattice. 
It leads to an electric potential both inside and outside 
of the GaAs substrate that couples to the electrons in 
graphene. The Hamiltonian of such piezoelectric inter- 
action between surface acoustical phonons of the GaAs 
substrate and the masslcss Dirac fermions in graphene 
can be written as 

H™ h = e^(R, t) = JL £ ^ A e< ,a 6q + cx , (1) 

where the electric potential ip(R,t) in the plane of the 
graphene sheet is given by the solution of the Poisson 
equation as p(R,i) oc q^e 1 ^'™^ e~ qd [5J]. Here 
q and u q — vpAq are the phonon momentum and en- 
ergy with the velocity of surface Rayleigh waves vpa ~ 
0.9v b 2.7 x 10 3 m/s, v b is the bulk sound veloc- 
ity in GaAs and q X}V = q x .y/q- In Eq. ([I]) 6 q denotes 
the amplitude of the phonon field and the piezoelectric 
electron- phonon interaction vertex is defined as = 

c pa (Qx<ly) 2 e~ 2qd h 2 vpA/ (poTpa) where we use the nom- 
inal time fpA ~ 8 ps, introduced in Ref. [35] for carrier 
scattering from bulk piezoelectric acoustical phonons. 
Here the characteristic wave vector po = 2.5 x 10 6 cm -1 
is related to the optical phonon energy in GaAs and the 
numerical factor cpa ~ 4.9 is determined by the elas- 
tic properties of GaAs [34]. The strongest electron-PA 
phonon interaction takes place for surface phonons prop- 
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agating along the diagonal direction with q x w q y so that 
in (|lj) we can approximate {q x %) ~ 1/4- Taking also 
e -2qd ^ £ or typical values of d ~ 5 A [301 135) one can 
see that after these well justified simplifications, the PA 
phonon vertex ^ A becomes independent of the phonon 
momentum q. This differs from the linear wave vector 
dependence of the DA interaction vertex and results in a 
new, qualitatively different, contribution to the mobility. 

The Hamiltonian of electron-acoustical phonon inter- 
action due to the deformation potential in graphene can 
be written as 



Hda = 



1 



V -y DA e 



i R 6 q + c.c. 



(2) 
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with the DA interaction vertex defined as |7~ 
H 2 qvDA/ (pqTda)- Here we introduce the nominal scat- 
tering time 1/tda = ^ 2 Po/ (^P v da) f° r DA phonon in- 
teraction and find fp>A ~ 0.8 ps using the following values 
for the graphene parameters [5]: S = 6.8 eV the coupling 
constant of the deformation potential, p = 7.6 x 1CP 7 kg 
vaT 2 the surface mass density of graphene, and vda — 
2.0 x 10 4 m/s the sound velocity in graphene, which is 
essentially larger from the sound velocity in GaAs. The 



2.0 x 10 5 cm- 1 /? 



direct comparison \"/q A \ / \lq A \ 
shows that both the PA and DA electron-phonon scat- 
tering mechanisms can be important in typical experi- 
mental situations and their actual contributions to the 
graphene mobility are determined by the carrier density, 
n, and temperature, T. 

The momentum relaxation rate of a test electron due to 
scattering from the potential of the s = PA, DA phonon 
field is 

ZI7^= E (1 - cos 9^)W^ x , k , 1 - ffi'fcfo ) 

where e>k = Xv F k is the electron energy in monolayer 
graphene with the chirality A and the Fermi velocity v F . 
The Fermi functions f{exk) are determined by the elec- 
tron Fermi energy e F and the lattice temperature T. The 
electron transition probability due to the emission (+) 
and absorption (— ) of phonons is 
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(4) 



with the Bose factors N (w sq ) representing the num- 
ber of phonons of the s,q mode. Making use of 
the electron wave functions in graphene V T (R) = 



Here # k k' = 6w — #k and the form factor J-A'A(#kk') = 
(1 + A A' cos #kk') /2 represents an overlap of the electron 
spinor wave functions. We consider electron-phonon scat- 
tering in doped graphene with e F much larger than the 
typical acoustical phonon energy. In such samples, only 
intra-chirality subband transitions are effective so we take 
A = A' = 1 and further omit this index. Then, we can 
represent the relaxation rate due to extrinsic PA and in- 
trinsic DA phonon scattering as 
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Q 2 +m{x) (6) 



where we introduce the functions 



f z± 

Gk(x) = ^2 dzz k r){a,z)^ ± {x,y) 
± Jo 



(7) 



with z = q/2k and z± = 1/(1 ± a). In 
Eq. ([6| for PA phonons we take m = and 
a = apA = vpa/vf, c = cpa', for DA phonons 
m = 1, a = a DA = v DA /v F , c = c DA = 2^2. 
For brevity, we define also the functions (x,y) — 
(N(y) + 1/2 ± 1/2) (1 - f{x T y)) I (1 - f(x)) where 
x = (e k - e F )/T and y = w sq /T = z(k/k F )(Tg A ' DA /T). 



The characteristic Bloch-Griineisen (BG) temperatures 
Tgg' DA = 2hvpA.DAk F are defined, respectively, for 
PA and DA phonon scattering. In Eq. Q the function 

r](a,z) = \Jl — (a ± (1 — a 2 ) z) 2 /(l T 2az) appears due 
to the chiral nature of Dirac fermions and restricts 
backscattering processes in graphene. As far as the 
Fermi velocity in graphene v F rs 1.15 x 10 6 m/s is 
much larger than the sound velocity in GaAs and 
in graphene, in our analytical calculations we take 
rj(a, z) m 77(0, z) = \f\ — z 2 and z± « 1. 

First we analyze analytically the momentum relaxation 
of a test electron in the regimes where PA and DA phonon 
scattering are qualitatively different and calculate the 
extrinsic PA versus the intrinsic DA phonon contribu- 
tions to the relaxation rate. For a typical doping level 
in experiments, n — n x 10 12 cm" 2 with n ~ 1, we 
have e F ~ 1350 K and for energies and temperatures 
up to room temperature, the system of massless Dirac 
Fermions is a well defined Fermi gas where carrier scat- 
tering events due to PA and DA phonons are kinemati- 
cally quasielastic, i.e. \e — e F \ , T <C £ F - Therefore, we 
consider three typical temperature regions with bound- 
aries given in terms of two different BG temperatures, 

8, all three regions are 



T. 



PA, DA 



bg ■ Because T^/T, 
well defined. 



PA 
BG 



ik-R 



In the high temperature regime, T ^ T BG ' , PA 
and DA phonons with energies u; sq ~ T BG are important 
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/V2A where k is the polar angle of the in the relaxation processes and because w sq < T, such 

scattering events are called [35] statically quasielastic. 
Under this severe condition, we have y <C 1 so that the 
phonon Bose factors are large, N(y) « l/y 3> 1, and 
scattering is dominated by induced phonon emission and 



electron wave vector k, we obtain for the square modulus 
of the PA and DA matrix elements 



Jk',k=pq 
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Al^kk' 



(5) 
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absorption processes. Therefore, in Eq. Q to leading 
order in y, we can replace f(x±y) by /(x) while N(y) + 
\ ± \ by l/y. Then, g 2 + m (x) « bT /Tbg (with m = 
and 6 = 1/3 for PA and m = 1 and 6 = 7r/16 for DA 
phonons) and for the momentum relaxation rate in this 
regime, we obtain 



be 1 
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PA. DA i x 



27T Tpyl^A \ Po 



, \ m+1 m 



(8) 



Here the relaxation rate due to DA phonons reproduce 
the previous results from Refs. [5] and The linear T 
dependence of the momentum relaxation rate both due 
to PA and DA phonon scattering and its independence 
on the energy e — Ef of a test electron are distinctive 
features of statically quasiclastic electron-phonon scat- 
tering in the high T regime where the scale of varia- 
tion of Ti(e) is Tp, which is larger than T BG , and the 
Pauli exclusion principle does not play an essential role. 
In this regime the extrinsic PA and intrinsic DA con- 
tributions to the momentum relaxation rate differ by a 
factor Tj DA (ek)/T ] PA (ek) = 7 HPo/kF ~ y/0.5/n where 
= bpAC %A fDA « 0.5 is determined by the elastic 
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properties of graphene and GaAs. Thus, in this regime 
we find that independent of e — Ef and T, the extrin- 
sic PA phonon contribution dominates the intrinsic DA 
phonon contribution to the relaxation rate at densities 
smaller than rt = 5 x 10 11 cm~ 2 and this enhancement is 
a square-root function with n. 

In the low temperature regime, T < TjgQ , if addi- 
tionally the energy of a test electron is small, \e — ep\ <C 
Tqq' DA , electronic transitions are dominated by small 
angle scattering events. The typical phonon momenta, 
q <C hp, correspond to the z — > limit and one can ex- 
tend the integration over z in Eq. Q up to infinity. Re- 
placing also 1 — z 2 by 1 and taking into account that z rj 
y (T /Tbg), we have Q k {x) (T /T B a) k+1 Fk{x) where 
we define the function Fk{x) = J2± Io° dyy kl $ ,± (x, y). 
Hence, the momentum relaxation rate in the low T 
regime is given as 
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Here one should distinguish two subregimes. For thermal 
electrons \s — £f\ < T, electronic transitions with typical 
phonon momenta hq ~ T/vp^ p>A dominate. The func- 
tion Tk{x) weekly depends on x for x < 1 and can be 
replaced by J-fc(O). Hence in this BG regime the mo- 
mentum relaxation rate Q exhibits T 3 and T 4 depen- 
dences, respectively, for PA and DA phonon scattering. 
For hot electrons in the opposite T40 limit, electronic 
transitions are dominated with spontaneous emission of 
phonons. We have Fk{x) « |x| fc /(l + k) for \x\ ^> 1 and 
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FIG. 1. The electron momentum relaxation rate 
l/T^ A ' DA (ek) as a function of temperature for different values 
of the electron energy and density. The dashed and dotted 
lines correspond to the extrinsic PA and intrinsic DA phonon 
scattering mechanisms. The solid curve is the total rate. Four 
different situations are presented: (top, left) Ek — £f = 3 K 
and n = 0.5 x 10 12 cm -2 , (top, right) Ek — £f = 3 K and 
n = 10 x 10 12 cm" 2 , (bottom, left) e k - ef = 15 K and 
n — 1.0 x 10 12 cm -2 , and (bottom, right) Ek — Ef = 30 K and 



the relaxation rate ^ does not depend on T and is pro- 
portional to ( \e — ef 

\l T BG DA ) m -The direct compar- 
ison of the extrinsic PA and intrinsic DA contributions 
to the momentum relaxation rate in the low T regime 
gives r 1 DA (e k )/r 1 PA (£ k ) = lL (p /k F ) (Tpa/T) for ther- 
mal electrons where the factor 7r = a ? AC ? A ^^!}j- DA ~ 
41. For hot electrons 7l is even larger. Thus, in the low T 
regime, the intrinsic DA contribution is negligible while 
the extrinsic PA contribution leads not only to strong en- 
hancement of the momentum relaxation rate but changes 
its dependence on the lattice temperature from T 4 to T 3 . 

In the region of intermediate temperatures, T^q <C 
T <C TgQ, small angle scattering events dominate the 
DA phonon relaxation with the rate given by Eq. Q 
while the PA phonon relaxation is governed by quasielas- 
tic large angle scattering events with the rate Q. There- 
fore, the momentum relaxation of thermal electrons due 
to scattering from intrinsic DA phonons is still suppressed 
in comparison with that from extrinsic PA phonons by a 
factor T? A tsK)/T[ A {e±) = 7l (p /k F ) (T^/T) 3 where 



bpAC pA TQA 



As seen this suppression depends 



^3(0)c 2 DA Tp A ' 

strongly on T, however, due to the small numerical pre- 
factor, 7/ w 4 x 10~ 3 (2.5 x 10~ 2 for hot electrons), it is 
large only near the lower edge, T > T^q, of this interme- 
diate temperature region with the relaxation rate linear 
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FIG. 2. The exponents a(n, T) (top panels) and f3(n, T) (bot- 
tom panels) describing the mobility behavior /i 3 oc T -a ("> T ) 
and /i s oc n /3 ( n ' T \ respectively, as a function of T and n. The 
panels on the left and right hand sides depict, respectively, 
the PA and DA phonon limited mobilities \ipa and fiDA- 



in T . Towards the upper edge, T < Tgg, the DA phonon 
contribution to the momentum relaxation rate increases 
and at high densities the PA contribution prevails with 
the T 4 behavior of the relaxation rate. 

Within the Boltzmann transport theory, the electron 
momentum relaxation time averaged over its energy e, 
rf(n,T) = JdeD(e)Tl(e) (-df(e)/de) (D{e) is the elec- 
tron density of states) , determines the carrier mobility in 
graphene samples as [T7] 



fJ,pA,DA(n,T) = 



ev F PA, DA 

hk F Tl 



(n,T) 



(10) 



As seen from Eqs. {[sj) and ^ in the high T regime as well 
as for thermal electrons in the low T BG regime, the mo- 
mentum relaxation rate of a test electron is independent 
of its energy. Therefore, we find that in the high T regime 
the mobility shows the same, \i oc T , temperature 
dependence for both the PA and DA scattering mech- 
anisms but different density dependences, /i cx l/y/n and 
fi oc 1/n, respectively, for PA and DA phonon scattering. 
In the low T BG regime the mobility exhibits the same 
density dependence, \i oc y/n, but different temperature 
dependences, fi oc T~ 3 and /i oc T~ 4 , respectively, for PA 
and DA phonon scattering. 

In Figs. [T]3] we present our numerical calculations 
based on Eqs. ^ and Jl0| ). In Fig. [I] we study the 
electron momentum relaxation scattered due to the com- 
bined action of extrinsic PA and intrinsic DA phonons. 
The relaxation rate versus T is plotted for different val- 
ues of e — £f and n. It is seen in all panels that at low 
temperatures the relaxation rate is almost independent 
of T, which corresponds to the behavior of hot electrons 
with e — £f ^> T. In the small angle scattering sub- 
regime, T< e-£f < Tqq' DA , the intrinsic DA scatter- 
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FIG. 3. (left) The mobility ratio ^pa/^da and (right) the 
total mobility \x = iipa[ida/([ipa + ft da) versus tempera- 
ture. The solid, dashed, and dotted curves correspond to the 
electron density n = 10, 5 and 2 x 10 11 cm -2 . In the right 
panel, the curves in each set are calculated in three differ- 
ent approximations for the PA phonon potential. The lower 
curves refer to the (q x q x ) 2 ~ 1/4 and d = approximation, 
the middle curves take into account the angle-dependence of 
the PA vertex with d = 0, the upper curves include both the 
effect of anisotropy and of the finite distance with d = 5 A. 



ing is strongly suppressed with respect to the extrinsic 
PA scattering (c/. the (top left) and (top right) panels) 
and this is in agreement with the above discussion that 

-, / PA.DA ( ir r PA,DA\ 3+m ■ ■ , 

1/Tj oc Ie — ep/l BG I m this subregime. 

With an increase of e — Sp the dependence weakens and 
in the large angle scattering subregime, e—ep > T^q' DA , 
the PA and DA phonon contributions to the momentum 
relaxation of hot electrons become of the same order and 
the DA phonon contribution can even dominate at low 
densities (cf. the (bottom right) panel). In the high T 
regime, as seen in all panels of Fig. [T] the momentum re- 
laxation rate exhibits the same linear T dependence, in 
agreement with Eq. ([8]). The PA and DA phonon mech- 
anisms make equal contributions to the relaxation rate 
at carrier densities n = 5.0 x 10 11 cm~ 2 (cf. (top left) 
panel). In this regime the relaxation rate is indepen- 
dent of e — £p and the relative PA and DA contribution 
can be tuned by changing solely the density n. At inter- 
mediate temperatures the momentum relaxation shows a 
crossover from the T independent to the linear T regime. 
The total rate depending on n and T can be dominated 
either by PA phonon scattering with the T 3 behavior for 
smaller values of n or by DA phonon scattering with the 
stronger T 4 dependence for higher values of n (cf. (top 
right) panel). 

In Fig.[2]we plot the exponents a(n, T) and /3(n, T), i.e 
1^pa,da oc T~ a ("' T ) and ^pa.da oc n /3 '™' T ', respectively, 
as a function of T and n. It is seen that the temperature 
and density behavior of the mobility calculated numeri- 
cally for the PA and DA scattering mechanisms is consis- 
tent with the above analytical findings. For any value of 
n, the T^ 1 and T~ 4 dependence of the mobility, respec- 
tively for PA and DA phonon scattering, occurs in a wider 
temperature range than the T~ x and T~ 3 dependence for 
DA and PA scattering. Such a reversed behavior stems 
from the substantial difference of the sound velocity in 
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GaAs and graphene and can be used to distinguish be- 
tween the extrinsic PA and extrinsic DA mechanisms of 
the momentum relaxation. As seen from the bottom pan- 
els in Fig. [2] at low T the exponent (3(n, T) ks 1/2 for both 
the PA and DA mechanisms. With an increase of T the 
PA contribution to the mobility shows a crossover to the 
behavior with f}(n, T) « —1/2 while the DA contribution 
exhibits a crossover to the behavior with j3(n,T) ps —1. 

In Fig. [3] we plot the temperature dependence of the 
relative and combined contributions to the mobility made 
by the PA to DA phonon scattering for different values 
of the electron density. It is seen from the left panel that 
at densities n < 5 x 10 11 cm~ 2 , the ratio ^pa/uda < 1 
so that PA phonon scattering is the dominant mecha- 
nism, limiting the mobility both at low and high temper- 
atures. The right panel shows the total mobility versus 
temperature, which we obtain applying the Matthiessen 
rule, l//x = 1/hpa + VM-DA- Here we calculate the PA 
contribution to the mobility including both the angle- 
dependence of the piezoelectric potential and the finite 
distance d between the graphene sheet and the substrate. 
From the comparison of the upper thick curves with the 
lower two curves in each set of lines, it is seen explicitly 
that the effect of the finite distance and the anisotropy 
of the PA potential is weak. 

In conclusion, the piezoelectric potential of acoustical 
phonons propagating on the surface of the GaAs sub- 
strate is an important factor in limiting the mobility of 
Dirac fermions. At low densities PA phonon scatter- 
ing is the dominant momentum relaxation mechanism 
in graphene. At high temperatures it changes qualita- 
tively the density dependence of the mobility while in 
the Bloch-Gruneisen regime the power law dependence 
on temperature. 
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